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The asymmetric nature of the fluid film stiffness and damping properties In 
rotors supported on fluid film bearings causes a forward or a backward whirl depend- 
ing on the bearing parameters and the speed of the rotor. A rotor was designed to 
exhibit backward synchronous whirl. The rotor-bearing system exhibited split crit- 
ical, and a backward whirl was observed between the split crltlcals. The orbital 
diagrams show the whirl pattern. 


INTRODUCTION 

The design of a rotor system must consider several aspects such as critical 
speeds, peak unbalance response, regions of change of whirl directions, and Insta- 
bility. In general, large rotor systems In continuous operation are supported on 
hydrodynamic bearings. These hydrodynamic bearings exhibit asymmetric cross-coupled 
stiffness and damping properties that vary with the speed of operation. Such a 
property Influences the dynamic behavior of the rotors significantly. 

The dynamic behavior of such rotors can be predicted by using transfer matrix 
methods, finite elements, modal analysis, etc. [1-6]. Rao [7] and Rao et al . [8] 
used a simple analytical technique to predict the dynamic behavior of such rotors. 
They studied a single-disk rotor on fluid film bearings and observed that for spe- 
cific rotor-bearing parameter combinations the system may exhibit two distinct peaks 
In the response but sometimes It may show only one peak In the response. Kollmann 
and Gllenlcke [9] have shown experimentally the existence of the split crltlcals In 
a simple rotor supported on fluid films bearings. Kellenberger [10] derived equa- 
tions for double frequency accelerations In turbogenerator rotors resulting from 
anisotropy In the plain cylindrical bearings and showed the occurrence of backward 
whirl between the crltlcals. These studies predict a forward whirl before the first 
critical and after the second critical, whereas the rotor executes a backward whirl 
between the two crltlcals. Also, the simple rotor with only a single peak In Its 
response Is known to whirl In the forward direction at all speeds. 

Several experimental Investigations are reported In the literature regarding 
the dynamic behavior of different types of rotor systems supported on hydrodynamic 
bearings. A few notable works are by Yamamoto [11], Hull [12], and Lund and 
Orcutt [13]. 

Most of the referred work on practical rotors supported on hydrodynamic bear- 
ings does not satisfy conditions for a clear backward whirl and hence the phenomenon 
of backward whirl has so far not been observed experimentally. The occurrence of 
backward whirl Is also not desirable In practice. In the present work conditions 
are derived for backward whirl considering bearing damping. Further results of the 
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experiments carried out on a laboratory rotor model designed so as to exhibit 
backward whirl are reported. 
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disk eccentricity 
shaft damping 

fluid film damping coefficients 

forward and backward component, respectively, of I** 1 modal force 
vector 

overall exciting force vector 
overall stiffness matrix 
shaft stiffness 

fluid film stiffness coefficients 
overall mass matrix 

forward and backward component of l^ modal displacement vector 
displacement vector at disk location 
displacement vector at bearing location 
maximum unbalance response of rotor 

displacement of rotor at disk location In Z and Y directions 

displacement of rotor at bearing location In Z and Y directions 

ith modal displacement vector 

I** 1 complex eigenvalue 

right eigenvector of system 

left eigenvector of system 


ANALYSIS 


A schematic of the rotor Is shown In Fig. 1. Equations of motion of the rotor 
are 


m 



(z + a cos wt) + k(z - Zg) + c(z - Zg) = 0 


(la) 


mz - mz = 0 


(lb) 
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my - my 


= 0 


m ^~~2 (y + a sin wt) * k(y - y Q ) + c(y - y Q ) = 0 
dt 

and the constraint equations are 


(1c) 

( 2 ) 


k( z - z 0 ) + c(z - Zq) = 2k zz • Zg + 2c zz • Zg + 2k zy • y Q + 2c zy • y Q 

k(y - y 0 ) + c(y - y Q ) = 2k yy • y Q + 2c yy • y Q + 2k yz • z Q * 2c yz • z Q 
Equations (1) to (4) can be written In the form 

[M] {5} + [K] {q} = (F) 

where 
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2 2 

{F } = [q, q, maw cos wt, maw sin wt, 0, 0] 

The eigenvalues and eigenvectors of the system are obtained by solving the 
homogeneous form of equations (6) as shown below: 
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( 6 ) 


[M] { Q( t ) } + [K] {Q( t) } = 0 
Assuming a solution of the form 

{Q( t ) } = {<p} exp (Xt) (7) 

where {<p} represent the system eigenvectors, and substituting equation (7) In 
equation (8) the eigenvalue problem becomes 

X[M]{<p} + [K] {<p} = 0 (8) 

Because of the presence of asymmetric cross-coupled stiffness and damping coeffi- 
cients in the bearings, the matrices [M] and [K] are nonsymmetrlc, resulting In 
a non-self-adjoint system. Hence a conventional normal-mode method Is not possible, 
and It Is essential to consider that the blorthogonallty property of the modes of 
the original system are those of the transposed system to the uncoupled equations 
of motion. THe left eigenvectors are obtained by transposing matrices [M ] and 
[K] In equation (6). The eigenvalues of the original and transposed systems are 
Identical, but the eigenvectors are different. 

The solution of equation (5) Is assumed In the form 

{Q(t)} = [<p]{n(t) } (9) 

where [<p] contains the eigenvectors of the original system (which are called 
right eigenvectors). Introducing equation (9) In equation (5) and premultiplying 
the result by [<p*] T , which Is the left eigenvector of the system, lead to the 
following uncoupled equations: 

[y]tn(t) } * M{n(t)} = {a} (10) 

where 

[y] = [<p*] T [M][<p] 

Cic] = [<p*] T [K][<p] 
and 

{«} = [<p*] t {F) 

Equation (10) results In uncoupled equation of the form 

uinj(t) + ic^ni(t) = <^( t) (11) 

This Is solved by assuming steady-state solution of the form 

ni(t) = Ni exp (jwt) + exp (-jut) (12) 

and 

<J 1 ( t) = Ei exp (jut) + E, exp (-jut) 
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Substituting equation (12) In equation (11) results In 
(*i + jwvi) I'M exp (jot) + (k^ - joy^) exp (-jot) 

= exp (jot) + Ei exp (-jot) (13) 

Equating the coefficients of the forward and backward rotation terms, we obtain 


N 


1 



E 


1 


+ j«y^ 


(14) 


^1 = k^ - joy^ 

Hence the displacements are determined from equations (10) to (14) and then the 
unbalance response of the rotor Is obtained nondlmenslonally as 


where a Is the disk eccentricity. The complex eigenvalues are obtained from 
equation (6). 


Depending on the bearing parameters the rotor will have either split crltlcals 
or just a single peak In the fundamental critical speed region, corresponding to 
synchronous whirl. J[n the case of the split crltlcals in the synchronous whirl, the 
backward component N^ Is larger than the forward component Ni between the 
split crltlcals. 


EXPERIMENTAL RESULTS 

The details of the test rotor are given In table 1. It consists of a circular 
shaft with a circular disc at the center and supported on two Identical hydrodynamic 
bearings at the ends. These bearings are mounted on cast iron pedestals at the two 
ends, and In turn these pedestals are rigidly fastened to the support, which Is made 
of steel angles. The pedestals were Impact excited both In the horizontal and ver- 
tical directions, and the resulting acoustical response was measured close to the 
pedestals. A frequency analysis of the measured sound showed that the first peak 
occurred at a much higher frequency than the critical speeds of the rotor, Indicat- 
ing that the pedestals were rigid. The bearings are supplied with oil through a 
gravity feed. The unbalance response of the shaft Is measured In both the x and y 
directions by two proximity pickups. The signals from the pickups are fed to a 
twin-channel FFT analyzer, and the orbit diagrams are obtained with a x-y plotter. 
The direction of plotter pen motion Indicated the direction of rotor whirl. 

In designing the rotor such that it exhibits a backward whirl, the bearing 
coefficients to be used In equations (3) and (4) were taken from the results of 
Lund [14]. These results can also be found In [15]. Lund's results agree well 
with the experimentally determined bearing coefficients by Gllenlcke [9]. 
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The whirl orbits obtained experimentally are shown normalized with respect to 
the disc eccentricity In Figs. 2 to 4 for different rotor speeds. Orbital diagrams 
are shown for a bearing clearance of 1.8796x10“* m, since a backward whirl could 
be identified only for this case. The two critical speeds are 2200 and 2600 rpm. 
The rotation of the rotor Is In the counterclockwise direction and hence Figs. 2 
and 4 show forward whirls and Fig. 3 shows backward whirl. The speed corresponding 
to the backward whirl In Fig. 3 Is 2500 rpm, which falls between the two crltlcals 
at 2200 and 2600 rpm. Since the backward whirl between the crltlcals Is of Inter- 
est, a photograph of backward whirl motion of the rotor at 2500 rpm was taken from 
the FFT analyzer screen. Theory also predicts these whirl directions, and hence 
there Is a qualitative agreement between theory and experiments. 
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TABLE 1: DETAILS OF TEST ROTOR 


Shaft Diameter 
Shaft Length 
Weight of Disk 
Shaft Stiffness 
Bearing Diameter 
Bearing L/D Ratio 
Oil Viscosity at 25.5°C 
Unbalance of Rotor 
Bearing Clearance 


0.0222 m 
0.5080 m 
89 N 

8.9 x 10 b N/m 
0.0254 m 

1 c 

0.96 x 10“° Pascal sec. 

1.084 x 1G“ 4 kg.m 

1.8796 x 10~ 4 m (0.0074 in) 
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